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Abstract 

We study the problem of the phonon-induced electron-electron interaction in a solid. Start- 
ing with a Hamiltonian that contains an electron-phonon interaction, we perform a similarity 
renormalization transformation to calculate an effective Hamiltonian. Using this transformation 
singularities due to degeneracies are avoided explicitely. The effective interactions are calculated 
to second order in the electron-phonon coupling. It is shown that the effective interaction be- 
tween two electrons forming a Cooper pair is attractive in the whole parameter space. For a 
simple Einstein model we calculate the renormalization of the electronic energies and the critical 
temperature of superconductivity. 
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1 Introduction 



BCS-theory [GJ is one of the most successful and most popular theories in condensed matter 
physics. It is based on a Hamiltonian with an attractive interaction between the electrons. 
Phonons are not explicitely present in this model. The phonon frequencies enter in the explicit 
form of the phonon-induced effective electron-electron interaction. Most of the properties of 
superconductors can be understood using BCS-theory. But one of the problems of BCS-theory 
is that the critical temperature T c of the superconductor cannot be calculated. The famous BCS- 
formula T c = 1.13G exp(— 1/N(€f)V), where is the Debye temperature and V is the strength of 
the interaction, has to be used to determine N(€f)V from T c . If one uses for N(€f)V values that 
have been obtained using standard perturbative treatments of the electron-phonon interaction, 
the calculated critical temperature is too large. Reliable results for T c can be obtained in the 
framework of the Eliashberg theory ||. For a review of this theory and a critical discussion of 
various formulas for the critical temperature we refer to || |H, see also §. 

In their discussion of BCS-theory, Allen et al Q argued that the BCS-equation yields a wrong 
result for T c because retardation effects in the interaction are neglected. This point of view is 
possible in a formulation of the theory in terms of Green's functions like in the Eliashberg-theory 
0] or equivalently in a field-theoretic formulation of a model containing electrons and phonons, 
where it is possible to integrate out the phonons. In both cases one obtains an effective phonon- 
induced electron-electron interaction that depends on time or frequency. But in the standard 
formulation of BCS-theory one starts with a Hamiltonian. The effective electron-electron in- 
teraction Bardeen et al had in mind was the phonon-induced interaction of Frohlich || or 
Bardeen and Pines Q. In these approaches the phonon-induced interaction is obtained by ap- 
plying a unitary transformation to eliminate the electron-phonon interaction in the Hamiltonian 
in lowest order. It is clear that in a Hamiltonian formulation of the theory, the interaction cannot 
be frequency dependent, since the Hamiltonian is a hermitian operator. From the viewpoint of 
a Hamiltonian formulation the problem of BCS-theory with an interaction of the Frohlich type 
is that the interaction is singular and does not contain the correct energy scale. If one performs 
a single unitary transformation to eliminate the electron-phonon interaction, one attempts to 
treat all energy scales in the problem at once. This usually fails, even in perturbation theory. 
Instead one should break the problem into pieces, dealing with each energy scale in sequence. 
This is the usual approach of renormalization theory. 

About two years ago Glazek and Wilson proposed a new renormalization scheme for 
Hamiltonians, called similarity renormalization. The main idea of this approach is to perform 
a continuous unitary transformation that yields a band-diagonal effective Hamiltonian. In the 
effective Hamiltonian the dependence on the original ultra-violet cutoff is removed, but it will 
contain additional interactions. In the present paper we apply this renormalization scheme to 
the classical problem of interacting electrons and phonons, modelled by a Hamiltonian 

H = ^2e k : c\c k : + ^u q : b\b q : +^2(g q c\c k+q b\ + g* q c\ +q c k b q ). (1.1) 

k q k,q 

c\, and Cfc are the usual creation and annihilation operators for electrons, b„ and b q for phonons, 
and the colons denote normal ordering. We calculate the renormalized Hamiltonian to second 
order in the coupling constant g q of the electron-phonon interaction. It contains an additional 
effective electron-electron interaction, that is responsible for superconductivity. 

Recently, the construction of the effective phonon-induced electron-electron interaction has 
been studied by Lenz and Wegner S. They used flow equations for Hamiltonians, a method 



proposed by Wegner [10] to block-diagonalize a given Hamiltonian. They obtained an expression 
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for the effective interaction that differs from Frohlich's result |J, in fact it is less singular. If we 
write the induced electron-electron interaction in the form 

2 Yl Vk > k '>1 '■ C k+g C k'- q C k'Ck (1-2) 
k,k',q 

the coefficient for the interaction of two electrons forming a Cooper pair obtained by Frohlich is 

V k ,-k, q = \g q \ 2 j- 2uJ \ 2 o (1-3) 

whereas Lenz and Wegner obtained 

V k ,-k,q = ~ \9q\ 2 ~ ( 2 - (1-4) 

There are two remarkable differences between these results. First, the interaction obtained by 
Lenz and Wegner has no singularity, second it is attractive for all k and q. 

Recently Wegner's flow equations for Hamiltonians have been applied to the Anderson im- 
purity model [Tll| . For that model a similar problem occurs. It can be mapped onto the Kondo 
model using the well known Schrieffer- Wolff transformation |T^| . The Schrieffer- Wolff transfor- 
mation eliminates the hybridization between the impurity state and the electronic band states to 
first order. Using this transformation, the induced spin-spin interaction contains a singularity if 
the impurity orbital lies in the electronic band. Furthermore, the induced interaction does not 
have the correct energy scale built in, therefore the Kondo temperature can not be calculated 
correctly. The flow equations yield a less singular expressions which has a different sign in some 
regions of the parameter space. In the case of the Anderson impurity model it was possible to 
compare the results with renormalization group calculations, and it could be shown that the flow 
equations yield the correct behaviour. For instance, the energy scale which sets the scale for the 
Kondo temperature comes out correctly, whereas the Schrieffer-Wolff result does not have this 
property |Q. 

Unfortunately, a similar comparison with renormalization group calculations was not possible 
in the case of the electron-phonon coupling. The primary purpose of the present paper is to fill 
this gap. Using the above mentioned similarity renormalization scheme, we will calculate the 
effective electron-electron interaction. We will show that the interaction between two electrons 
forming a Cooper pair has the sign structure of the expression obtained by Lenz and Wegner, 
it is attractive for all k and q. This clearly supports the result by Lenz and Wegner ( |1.4[ ). The 
difficulties of the T c formula of BCS-theory can be traced back to the wrong effective electron- 
electron interaction. We will come back to this point in Sect. 4. 

For a reader who is familiar with Eliashberg theory, it is probably unclear how an approach 
using a Hamiltonian framework can lead to quantitative results without a fit parameter like in 
BCS-theory. Neither the interactions (|l.3| ) or (|l.4| ) nor the interaction obtained using similarity 
renormalization contain retardation effects. The reason that the procedure yields reliable results 
is that similarity renormalization yields the correct low energy scale of the problem. Different 
energy scales are automatically separated. How this happens will be made explicit in Sect. 2. A 
second point is that since the renormalization procedure consists of a continuous unitary transfor- 
mation applied to the Hamiltonian, one has to transform the observables as well. Retardation, 
i.e. the decay of single particle states and the broadening of peaks in spectral functions can 
be obtained within the present approach, if the transformation of the observables is taken into 
account. In a recent paper on time-dependent equilibrium correlation functions in dissipative 
quantum systems it has been demonstrated, how the transformation of the observables can be 
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done [13]. In that case Wegner's flow equation [10] have been used, but the transformation 
of the observables is very similar in both approaches. It was possible to obtain very accurate 
quantitative results in that case. In the present paper we only calculate static properties of the 
system, therefore we will not discuss this point further. 

The outline of this paper is as follows. In the next section we give a brief description of the 
similarity renormalization scheme. For details we refer to §. This method has been applied 



to quantum electrodynamics [14] and to quantum chromo dynamics [15, [l6| (and the references 
therein). Many details concerning the method are explained in these papers as well. In Sect. 
3 we apply the method to the electron-phonon coupling. We derive the renormalization group 
equations and calculate the renormalized Hamiltonian to second order in g q . Mainly for illus- 
trational purpose we derive some explicit results for the simple Einstein model for phonons in 
Sect. 4. First we calculate the renormalization of the electronic single particle energies. They 
show a typical logarithmic renormalization, which leads to the well known renormalization of 
the density of states at the Fermi surface. Next we calculate T c within the Einstein model using 
our effective interaction. We compare our results with those obtained using ( |1.4| ) or using the 
T c -formula of BCS-theory. Finally we give some conclusions including a comparison between the 
similarity renormalization scheme and Wegner's flow equations. 

2 Similarity renormalization 

This method yields renormalization equations for a given arbitrary Hamiltonian H\ as a function 
of the cutoff A. The goal of the method is to transform the initial Hamiltonian (with a large 
cutoff A) into an effective Hamiltonian that has no matrix elements providing energy jumps 
large compared to the small cutoff A. The initial Hamiltonian contains usually multiple energy 
scales and couplings between these energy scales. The renormalized Hamiltonian H\ has a band 
diagonal form. The different energy scales in the initial Hamiltonian are decoupled in H\. 
Let us write H\ in the form 

H\ = Hq\ + Hj\. (2.1) 

Hq\ is the free Hamiltonian, Hj\ contains interactions and counter terms. The eigenvalues of 
Hq\ are Ei\ > 0. We introduce cutoff functions Uij\ for the matrix elements. To assure a band 
diagonal structure of the renormalized Hamiltonian, we choose Uij\ = 1 if \Ei\ — Ej\\ is small 
compared to lambda and U{j\ = if \Ei\ — Ej\\ is large compared to A. The detailed form of 
Uij\ is not important. One possible choice, which has been used in the treatment of QED on the 
light front p3 is 



u ijX = e(X-\E iX -E jX \). (2.2) 
In other situations one has to introduce a smooth cutoff function. For a discussion of this point 



we refer to [||, 15 1. Furthermore we introduce 

= 1 - Uijx (2.3) 
The renormalization of H\ can be written down as an infinitesimal unitary transformation 

lt = [m > Hx] (2 - 4) 
with a generator r]\. rj\ has to be chosen so that the matrix elements of H\ obey 

Hij\ = Uij\Qij\. (2.5) 
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This assures that H\ is band diagonal. For the matrix elements Qijx we obtain 

^j^QijX + = VijxiEiX ~ Ejx) + [vx, H IX ]ij. (2.6) 

One defines 

Gij\ = f [r]\,Hix]ij c[\~Q i: > x 

and chooses 

r)ijx{Eix — Ejx) = —rijxGijx- (2.8) 

This yields the final equations for the matrix elements of the generator of the infinitesimal unitary 
transformation 

r ijX ( , „ i dUijX HijX 



and for the renormalization of the Hamiltonian 

- 5r = ^A,HiA] y +r yA - 5r — • (2-10) 

Since r^^ = if — ;$ A, the denominator £^ — Ej\ in ( |2.9| ) is bounded from below. 
These equations are exact, but it is clear that one is not able to solve them explicitely. Suitable 
approximations are necessary. Since small energy denominators are avoided in this approach one 
possibility is a systematic expansion in a coupling constant. 

Before we proceed to the application of the similarity renormalization, let us take a closer 



look at the equations (|2.9| ) and ( 2.10 ). The main point is that due to ( |2.q ) the renormalized 
Hamiltonian contains only couplings between states with an energy difference less than A. Start- 
ing with an initial cutoff A, couplings between very different energy scales are eliminated first, 
couplings between states belonging to the same energy scale are not eliminated or are eliminated 
later for smaller values of A. As a consequence, different energy scales are seperated during the 
renormalization process. 



3 Application to the electron— phonon coupling 

The Hamiltonian of the model is 

Hx = Hqx + Hix (3.1) 
with 

#0A = ^2 £fc,A : c\c k : + ^ w 9iA : b\b q :, (3.2) 

k q 

H ix = J2- 9k^xc\c k+q b q + g* k , q> x c l+q c k b q : +°(f 2 )- (3-3) 

k,q 

ct and Cfc are the electron creation and annihilation operators. The index k is a shorthand 
notation for (k, a) and — k denotes (— k, a). b q and b q are the creation and annihilation operators 
for phonons. Let us assume u qt x = ^>-q,x, £fc,A = £-fc,A- For a large cutoff A we let gk,q,A = 9q,A = 
9- q ,A independent of k. The terms 0(g 2 ) in Hjx contain induced interactions. For simplicity we 
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do not include a Coulomb repulsion in our model Hamiltonian. We have introduced a normal 
ordering for the fermions and for the bosons, 

c\c k := c\c k - n k , (3.4) 

: b\b q := b\b q - n q , (3.5) 

where n k is the occupation number for electrons in the state k and n q is the occupation number 
for the bosons, respectively. The Fermi energy is set to zero. 

The model contains several energy scales. The largest energy scale is given by the band width 
of the electronic single particle states, which is usually a few eV. Typical phonon energies are 
given by the Debye frequency, they are about two orders of magnitude smaller. Due to the phonon 
induced interaction between the electrons, Cooper pairs are formed. Their energy scale is given 
by the critical temperature of superconductivity, which is again about two orders of magnitude 
smaller than typical phonon energies. One has to resolve an energy scale in a Hamilonian that is 
dominated by energies being four orders of magnitude larger. Such a situation often occurs when 
a given problem has a marginal relevant operator. In our case this is the attractive phonon- 
induced electron-electron interaction. When this interaction is calculated within a perturbative 
scheme like the one used by Frohlich, all energy scales are treated at once. This usually fails. 

We now investigate this model using (2.9, |2,10 ). All quantities will be calculated as series 



expansions in g. Due to (|2.9| ) the generator of the infinitesimal unitary transformation can be 
written as 



^X = Yl '■ ^k,q,Xclc k+q b\ - Vk !q ,X C k+q C k b q : + (# 2 ) ( 3 - 6 ) 
k.q 



with r) kq \ = 0(g). In order to have a definite choice for u ks ^\, we let 

Uk,q,X = U(\e k ,\ ~ £fc+g,A + Uq,\\ /A) (3.7) 

with a smooth cutoff function u(x) that drops fast from 1 to zero in the vicinity of x = 1. Then 
(P) yields 



r k , q ,x d\nu k ^ x ,nf^ a s^ 

Vk, q ,x = — ; — -77 9k, q ,x + 0(g ). (3.8) 

Cfc.A — £k+q,X + ^q,X uA 

Let us now calculate the effective Hamiltonian. The first term in ( [2.10D is of order g 2 . This term 
contains the renormalization of the single particle energies e k) \ and ujq t \ and some new, induced 
interaction in Hj\. The second term in ( 2. 10| ) is of order g and contributes to the renormalization 



flow of the electron-phonon coupling. The renormalization of this coupling is determined by the 
equation 

d9k,g,x d\nu k ^ x . 3 , f . 

= r kt g, x -j^9k^x + 0(g ). (3.9) 

It can be solved using the ansatz 

9k, q ,x = g q ,K-^^- ( 3 - 10 ) 

with e k>qj x = e(|e fc ,A - tk+ q ,x + u q ,x\ /A). e(x) obeys = r^H. Using r = 1 - u it can be 
written in the form 

e(x) = u(x) exp(r(x)). (3-H) 
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e{x) has properties similar to u{x): e(x) = 1 if u{x) = 1, e(x) = if u(x) = 0, e(x) decays 
monotonously. Let us assume that the cutoff A is sufficiently large, so that, e^^A = 1 holds for 
all k, q. This yields g k , q ,\ = g g ,Ae k) q,x- 

Let us now calculate the terms of order g 2 . To do this we have to calculate the commutator 
[r)x,Hix]- It is given by 

k,q k',q' 

+ h.c. +0(g 3 ) 

= ^2Y1 r lk,q,\gk',q',\(:- 4: C k+q+q' b l b l> : S k',k+q- ■ A-q' C k+q h \ h \' : S k',k-q> 
k,q k' ,q' 

+{n k - n k+q ) : b\bl_ q : 5 q - q >6 k r )k+q ) 

+ ^ T lk,q,X9k',q',x( : c k c k+q-q'b\b q > I 5 k ' jk +q-q'- ■ c k+q > c k+qb\b q ' \ 5 k > >k 
k, q k', q ' 

+(n k - n k+q ) : b\b q : 6 q>q t6 k > tk 

~ '■ C k C k' + q C k' c k+ q ■ 5q,q' 

-(1 - n k + n q ) : c\ +q c k+q : 5 qs >5 k ^ k + {n k+q + n q ) : c\c k : 5 q) q>5 k > jk 

+{(n k - n k+q )n q - (1 - n k )n k+q )8 m >5 k > )k ) + h.c. + 0(g i ) (3.12) 

This expression together with ( |2.10| ) yields the renormalization of u qt \ and e^. For u qt \ we 
obtain 

~ ^2(Vk, q ,xg* k , q ,x + Vk, q ,\gk, q ,x)(n k - n k+q ) (3.13) 

k 

Using (|3.8[) and the definition of e k ^x yields 

_ gk, q ,x d\ne Kq ^ 

Vk, q ,x : "7T (A i4 J 

Cfc.A — £k+ q ,X + ^q,X «^ 

and therefore 

du q ,x _ n k+q -n k dh\e\ qX 

dX e fc ,A - efc+ g ,A + u q> \ dX ' q ' 

n k+q -n k a\g k , q ,x\ (3 15) 

^ £fc,A — efc+ g ,A + ^> q ,x dX 

We integrate this flow equation from a small cutoff A to the initial large cutoff A. On the right 
hand side we perform an integration by parts. Then the remaining integral contains derivatives 
of the single particle energies and is therefore of fourth order in g kq \. Neglecting this term we 
obtain 

u q ,X = w g ,A + V — — V^±l j^l 2 (! _ e 2 ), (3.16) 

*-r £fc,A - (-k+ q ,x + u q ,x 

A similar calculation can be performed for e k \. The renormalization of e ky \ is determined by 

de k ,\ 
dX 



^{Vk^Xgl^X^k+q + Uq) + Vk,q,Xgk,q,x{n k +q + Tl q ) 



1 



-Vk+q-q,xg k+ q,-q,x(l ~ n k+q + n q ) - r)l +q _ q>x g k+q ^ q> x(l - n k+q + n q )) 

i 1 2 1 2 

n k+q + n q d\g kA) \\ ^ y-^ 1 - n k+q + n q d \g k+q _ q ^\ (3 17) 

q £fc,A — e k+q,X + ^q,X dX ^ € k+q ^\ — € kt \ + UJg^x dX 
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and with the same assumptions as above we obtain 

. n k+q + n q 2 /-, 2 \ 

efc.A = £fc,A + > \9q,M (1 - e fc g A ) 

^ £fe,A - efc+g,A + W 9i a 



1 - + n g 

e k+q,X — £fc,A + ^g,A 



E 



|9,,a| 2 (1-4+,,-,.a) (3-18) 



The two equations (|3.1q ) and ( 3.18 ) can be used to determine the single particle energies self 
consistently. The result is correct up to second order in g. The main point is that due to the 
factors (l-e 2 k q X ) or (1 - 

fifc-l^—g^) small energy denominators are avoided explicitely. Without 
these factors and with the renormalized energies in the second and the third term on the right 
hand side replaced by the unrenormalized energies, the renormalization equations ( 3.16| ) and 



( 3, 18] ) are well known and can be found in various textbooks on solid state theory (see e.g. fll8|). 



In the next section we show that ( |3.18| ) yields the correct renormalization of the density of states 
at the Fermi surface, which is usually calculated within the framework of Eliashberg theory. 

[r]\,Hfx] contains further terms generating new couplings in the Hamiltonian Hj. It can be 
written as 

HlX = 9k,q,\c\c k+ qb\ + 9k,q,X C k + q C khq : 

k,q 

+ \ Yl Vkk '4> X '■ C k+q C k'-q C k' C k ■ 
k,k',q 

+ Y.( C l b l b -4 + C 1 b 4 b -^ 
1 

+ couplings between electrons and two bosons + 0(g 3 ). (3.19) 

Let us calculate Vk,k',q and c q to second order in g. The other couplings can be obtained similarly. 
Following ( 2.10| ) the renormalization equation of Vkk'q,x is 



fc ^' 9 ' A = -Uk,k',q,x(Vk+q,-q,Xgi'-q : X + Vk' -q,X9k+q-q,X + Vk' -q,q,X9*k,q,X + Vk,q,X9k' ,-q,q,x) 

+rk : k',q,X ^ Vk,k',q,X- (6.ZU) 

In the first term in the right hand side the induced electron-electron interaction is generated due 
to the elimination of the electron-phonon interaction. The second term eliminates the induced 
electron-electron interaction again. This elimination then yields higher interactions which are of 
fourth or higher order in the electron-phonon coupling. The problem is now that we expect to 
obtain an attractive induced interaction between the electrons. Such an attractive interaction is 
known to be a marginal relevant operator for a fermionic system. The elimination of such a term 
can cause difficulties. Within the framework of Wegner's flow equations, this has already been 



observed for the one-dimensional problem [ 10 1 . Therefore we modify the renormalization scheme 
at this point. We simply choose u^k' q x = 1- Then the second term in the renormalization 
equation for the induced interaction vanishes. Using the above expression for gk, q ,x and %,<j,a> 
this equations becomes 

dVk t k',q,\ , ,2 ( e k'-q,X dek+ q -q,X &k+q-q,X ^ e fc',-g,A 

\9qM ~ — ; — — T\ h 



dX ' V e £i+9,A — e k,X + ^q,X d\ Ek'X ~ ^k'-q,X + ^q,X dX 

_| e k,q,X de k i-g,q t X | &k'-q,q,X de k ^ t x \ (3 21) 

tk>-q,X — £fc'A + ^q,X dX €k,X ~ ^k+q,X + ^q,X dX J 
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Defining 

fk,k', q ,x = J dsek'-q, q , s , (3.22) 

we can write Vk t k',q,\ i n the form 



Vk,k',q,X — — \9q,\\ 2 ( ~ I 1" 



fk,k',q,X fk'—q,k+q,q,X 



£fc,A — £fc+g,A + ^q,X ^k'-q,X ~ £fc',A + ^g,A 
fk+q,k'-q,-q,X _^ fk',k,-q,X ] _j_ Q^4\ (3 23) 



e fc+<?,A — Cfc,A + ^g,A ^k',X ~ £k>-q,X + ^,j,A, 

The terms 0(g 4 ) arise again due to the fact that the derivative of the single particle energies is 



of order g 2 . Further corrections 0(g 4 ) occur since in ( 3.21 ) terms of 0(g 4 ) have been neglected 



( 3.23; ) is our main result. Note that fk,k',q,x > 0. Let us discuss two interesting cases: 



a) k' = k + q (the diagonal part of the interaction) 
We have fk,k+ q , q ,x = g(l - e t, q ,x)- This y ields 



\r u i 2 / 2w 9> A 

v k,k+q,q,X ~ \9q,A 



(efc+g,A - £fc,A) 2 - w£ A 



=2 p 2 



e fc+c?,A — £fc,A + ^q,X £fc,A ~ £fc+<J,A + ^q,X 



(3.24) 



The first term is the well known expression already obtained by Frohlich || . The singularity 
is cancelled by the two other terms. These terms occur because the electron-phonon 
coupling has not been eliminated completely. 

b) k' = —k (the interaction of two electrons forming a Cooper pair) 
Using fk,-k,q,x = f-k,k,- q ,x we obtain 

t t o I „ |2 / fk,-k,q,X fk+q,-k-q,-q,X \ /, ,r\ 

Vk-k,q,x = -2 \g q ,A\ ■ 1 ■ . (3.25) 

V £fc,A — £k+q,X + ^q,X e fc+q,A ~ £fc,A + ^q,X J 

In the case uj qj \ > {e^+q^x — £k,x\ both terms are negative, the interaction is thus at- 
tractive. If u> q) x < |efc+<j,A — £fc,A | the sign of the two terms is different. Let us con- 
sider the case €k+q,x — £fc,A > u q,x- Under this condition one has \ek+ q ,x — £k,x + ^<j,a| > 
|efe,A - £fc+ ? ,A + ^q,x\- For some lar § e cuto ff A, both, |e fc+(?i A - £fc,A + ^ q ,a\ and |e fc ,A - 
e k+q,A + w ?,a| are less than A. Lowering the cutoff, it happens that \ek+q,x — £k,x + ^q,A| 
becomes larger than A whereas \ek,x — £k+q,x + ^q,x\ is still smaller than A. This mean that 
fk+q-k-q-q,x is nonzero, whereas fk-k,q,x remains zero. In that case only the second term 
in ( |3.25 ) contributes. This term is negative, the interaction is thus attractive. Similarly, 



if £fc,A — efc+<j,A > ^q,x, the first term starts to contribute first and is negative. This shows 
that for sufficiently large A and Vk : -k, q ,A = one has Vk : -k, q ,x < 0. This property holds 
for general A as well if e{x) (or u(x)) falls of sufficiently fast in the interval around x = 1. 
To see this, we calculate the integral in the definition of fk,-k,q,\- The integration from A 
to A can be replaced by an integration from Ai to A, where Ai = A if ek, q ,xek+q-q,x > 0, 
otherwise Ai is the supremum of all s with c^ q gd^^q ^q g = 0. In a second step we use that 
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efc >g>s , and ek+ q -q )S ar e monotonic increasing functions of s. Then the second mean value 
theorem for integrals yields 

fk,-k,q,X = (i s dln ^,q,s = _ lnefcgA2 (326) 

with some A2 > Ai. This finally yields 

K = 2 ^ ( * + In^A (3 27) 

A2 V a; fc+q,-g,A 2 x k,q,X 2 J 

where Xk, q ,x = (efc,A — ^k+q,x+^q,x) / X. Whenever |lne(x)| / \x\ is a monotonously increasing 
function of \x\, the right hand side is negative. This condition can easily be satisfied. Let 
us write e(x) in the form e(x) = exp(— |x| Then the condition is satisfied if h(x) 

increases monotonously as a function of x. This means that the functions e(x) and u(x) 
have to decay sufficiently fast. 

To obtain an explicit expression for the interaction within a Cooper pair, one has to choose a 
specific form of e(x). If we assume that e(x) drops rapidly from 1 to 0, we can take for simplicity 
e(x) = 6(1 — x)6(l + x). This choice of a step function is useful for doing analytical calculations, 
but as mentioned above it can lead to pathologies in higher orders in g. With this choice of e(x) 
we obtain 

fk,-k,q,X = #(e/c,A ~ ^k+q,x)^(\^k,X ~ ^k+q,X + Vq,\\ ~ A). (3.28) 

This yields 

2 | 5 a |2 

Vk-k,q,X = — | ' j— ; 6(\€-k+q,\ — e k,x\ + ^q,X ~ A). (3.29) 

\£k+q,X — £fc,A| + ^>q,X 

The form fl3.29|) of the phonon-induced electron-electron interaction is similar to the result 
obtained by Lenz et al ||, see ( |1.4[ ), except for the cutoff function. It differs from the result by 
Frohlich (|Q|). The three expressions agree for \ek+q y x — £k,\\ = 0- 
The coupling c 9i a can be calculated likewise. We obtain 



diCq t \ 1 \ . 

^ = —2 U K X q,X) ^ n k{r]k,q,Xgk+q,-q,X + Vk+q-q,Xgk,q,X + r]k-q,Xgk-q,q,X + r]k-q,q,Xgk- q ,X) 
k 

+r(x q ,\) TT^^a (3-30) 



where 

Xq,X = (3.31) 
With a similar ansatz as above, 

c q ,x = e(xq,x)c q ,x (3-32) 
we have 

^rr- = -\^V{-r{Xq t x))^n k (ri KqtX gk+q-q,X + llk+q-q,Xgk, q ,X 



dX 

k 



+Vk,~q,Xgk~q,q,X + Vk-q,q,Xgk,-q,x) ■ (3.33) 
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c qi x can be calculated using the same approximations as above. We obtain 



L I 2 I I Jk,q,X 



fk,q,X . fk+q,—q,X 



£fc,A — £fc+9,A + ^>q,X £fc+(?,A ~ Cfc.A + ^>q,X 
_| fk,-g,X | fk~q,q,X \ ^ 



where 



fk, q ,x= I ds exp(-r(x qiS ))e k+g _ q>s dek > q > s . (3.35) 



A 



ds 



In the final Hamiltonian these couplings can be treated perturbatively. They yield a renor- 
malization of the phonon frequencies in fourth order in g g A- Similarly, the additional terms in 
the Hamiltonian Hj which contain couplings of an electron to two phonons yield an additional 
contribution to the electron-electron interaction of fourth order in g. 



4 Some explicit results for the Einstein model 

The above results are very general and can be applied to any type of electron-phonon interac- 
tion, even to any type of interaction between electrons and other bosonic degrees of freedom. In 
this section we apply our results to the simple Einstein model of phonons. Furthermore we will 
assume a constant density of states for the electrons near the Fermi surface, and we neglect renor- 
malization effects of the phonon energies. The motivation is to show the reader how well known 
results can be obtained using the above procedure. We will also show what are the main effects 
of the new phonon-induced electron-electron interaction ( 3.29Q on the critical temperature. 



In the Einstein model one takes phonons with u> q = ojq and g q ^\ = goV~ l l 2 . These assump- 
tions yield considerable simplifications in the above formulae. Let us first study the renormal- 
ization of the electronic energies e kj \. From ( |3.18| ) we obtain 

efc,A = £fc,A + glV" 1 y~] — ■ 0(\e k x - e fc '.A + w | - A) 

^ £fc,A - Cfe'.A + ^0 

V l - nk ' g(| e _ e + UJo \- A ), (4.1) 

^ efc'.A - £fc,A + 

at zero temperature, where the phonon occupation number n q = 0. As usual we replace the 
summation over k' by an integration over ey. Furthermore we assume that the density of states 
is constant, we denote it by N(0). The integral over ey extends over the whole band and we 
assume that the band width D is large compared to ujq and A. The renormalization of band 
energies ek,x with \e^\\ < wo is then given by 

, 2A7-/n\i ^0 - £fe,A + max(0, A + e fciA - ^o) (A n , 

ek,x = e k> A + g N(0)ln — — -. (4.2) 

wo + £fc,A + max(0, A - e k ,x - wo) 

For a cutoff A that is smaller than the phonon frequency luq and for e k <C wo this yields the well 
known expression (see e.g. Q) 

€k ' X = 1 + 2^(0)/^ (43) 

The last expression can be derived in the framework of Eliashberg theory, which is valid under 
the same conditions as used for the derivation of (|4.2j), namely ojq <C D. An important point in 
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this discussion is that the cutoff A should be somewhat less but of the order of ojq. With this 
choice properties of electrons close to the Fermi surface are described correctly. If the cutoff 
is larger, a part of the electron-phonon interaction that is relevant for the electronic behaviour 
near the Fermi surface has not been eliminated. A lower cutoff does not affect the single particle 
electronic properties near the Fermi surface. 

Let us now come to the discussion of the phonon-induced electron-electron interaction. The 
expression ( 3.29| ) for the interaction of two electrons forming a Cooper pair contains the factor 



$(l e fc+g,A — £k,\\ + ujq — \)- This factor arises due the fact that for a given A only a part of the 
electron-phonon interaction has been eliminated. If we choose A to be less than ujq, this factor is 
unity. This is similar to the behaviour of the single particle energies. This part of the interaction 
then no longer depends on the cutoff A. 

The energy gap and the critical temperature can be estimated if we treat the renormalized 
Hamiltonian as in BCS-theory. This leads to the usual BCS equation for the energy gap, where 
the electron-electron interaction ( 3.29| ) has to be taken, 



In the isotropic case we can replace the sum over q by an integral over et+q- Introducing 
A(efc) = A/% and using ( p. 29 ) we obtain 



A(e) = 2gl [ de' pit')- \ , A ^ tanh ( £ Je' 2 + A(e')A . (4.5) 

W °J K '\e-e'\+LU 2 v / e' 2 + A(e') 2 V2 V K 'J 



Here p(e') is the renormalized density of states since (4.4) has been derived from the renormalized 
Hamiltonian and therefore we have to use renormalized energies. It can be approximated by a 
constant renormalized density of states p(0) = iV(0)/(l + 2<7giV(0)/u;o). Using dimensionless 
quantities x = e/mo, x' = e'/ujQ, r = (/3u;o) _1 , A(x) = A(e)/uio we obtain 



A(x) 2 9o N (°)/^ f dx > A (^) tanh(Vx'2 + A(x')7T) , . 

y) 1 + 2g 2 N(0)/cj J \x-x'\ + l vV 2 + A(x') 2 

This equation can be compared to the usual BCS equation with a constant interaction in an 
energy interval given by ujq, 

■> , x , f ,tanh(vV 2 + A 2 It) , . 

1 = 2g 2 N(0)/u; / dx' \ 11 (4.7) 

J\x'\<i Vx' 2 + A 2 



The main difference between our result ( [4.6| ) and the BCS result is the renormalization of the 
density of states N(0). It is well known (see e.g. the discussion by Allen et al, ||, |5|) that this 
renormalization is important when the coupling 2g 2 l N(0)/uJo is not small. Another important 
point is the different way the cutoff is introduced. In BCS theory one has to introduce a cutoff 
by hand, and this is usually done by restricting the electronic energies to an energy interval of 
width 2coq around the Fermi surface. The main motivation for this is that an interaction of the 
Frohlich-type ( |1.3[) contains a singularity and becomes repulsive for larger energies. In our case 
we use a renormalization procedure which leads automatically to a cutoff ljq due to the factor 
(\x — x'\ + 1) in the denominator. 

(4.6) cannot be solved analytically, but it is easy to solve it numerically. In Fig. 1 we show the 



solution for A(x) for 2g 2 N(0)/uJo = 1/3 and 2gQN(0)/uio = 1.4 for different temperatures. One 
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clearly sees that the gap becomes small for large x and decreases with increasing temperature. 
The scaled curves in Fig. 1 show that the scaling relation 

A(„,T) = A(„,0)^|a (4.8) 

is satisfied within the numerical accuracy. We have checked this scaling law numerically for 
values of 2<?q./V(0)/u;o up to 2.5. The scaling law has been found by Scalapino et al [17] (see also 



||]) for weak coupling superconductors, whereas in the strong coupling case deviations occurred. 
But their calculations have been done for more complicated phonon spectra so that a direct 
comparison with our results is not possible. Scalapino || also mentions that even for strong 
coupling the reduced energy gap A(0, T)/A(0, 0) as a function of T/T c follows very closely the 
BCS curve for this quantity. We have plotted the reduced energy gap in Fig. 2, again for 
2gQiV(0)/wo = 1/3 and 2gQN(0)/u>o = 1.4. For other values of the coupling strength the curve is 
the same, it agrees with the BCS curve. The data shown in Fig. 2a show that for small values 
of the coupling the ratio 2A(0,0)/T C = 3.53 as predicted by BCS-theory. For 2g%N{0) = 1.4 
we obtain A(0,0) = 0.298w and T c = 0.161w - In that case the ratio 2A(0,0)/T C = 3.70. 
This tendency to larger values of 2A(0,0)/T C for strongly coupled systems in well known. The 
behaviour of A(u>,T) for other values of 2gQN(0)/u>o is similar. The main difference is the value 
of A(0), which shows a strong dependence on 2<7giV(0)/u;o as expected. 

An important point is that the energy gap A(0) is about an order of magnitude (for small 
couplings 2gQN(0)/ujQ ;$ 0.5 a factor of 5) smaller than the value one obtains using (fl.Tj). As a 
consequence the calculated critical temperature is smaller by roughly the same factor. We show 
this general feature in Fig. 3, where T c is plotted as a function of 2<?Q./V(0)/a;o together with 
the BCS-result. A comparison with the value of the critical temperature calculated with the 
interaction ( |1.4j ) from Lenz and Wegner || is not so easy, since in their paper the renormalization 
of the electronic single particle energies was not taken into account. But the flow equation for 
in their paper is similar to what we obtained, so that one should expect a similar renormalization 
of the single particle energies. Therefore we put this in by hand. The values for the gap calculated 
with the same assumptions as above show a behaviour that is very similar to our results. Only for 
larger values of 2gQN(0)/u!Q deviations occur. For 2gQN(0)/cuo si 1.2 the gap for the interaction 
(|1.4| ) lies below our curve. Furthermore it is interesting to compare our result with results for 
T c obtained in the framework of Eliashberg theory. One famous result is the McMillan-Dynes 
equation E(J], which has in our case the form 

T.-^P^t^y^V (4.9) 



1.2 K V 2g$N(0)/u, a I 

This curve has also been plotted in Fig. 3. The plot shows that for small and intermediate 
coupling the agreement is very good, whereas for strong coupling (2<?q N(0)/loo ^ 1.2) deviations 
occur. It is well known that in the strong coupling case the McMillan-Dynes equation is not 
valid. But in the region where the McMillan-Dynes equation is applicable, we obtain a good 
agreement. These results clearly show that a renormalization procedure for Hamiltonians like 
similarity renormalization or flow equations treats the different energy scales in the problem in 
a satisfactory way. 

5 Discussion of the results 

Starting from an initial Hamiltonian with an electron-phonon interaction, we calculated an effec- 
tive Hamiltonian to second order in the electron-phonon coupling. Due to the special structure 
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of the similarity renormalization scheme, the effective Hamiltonian has a band-diagonal struc- 
ture. It contains still a small part of the electron-phonon interaction, in addition an effective 
electron-electron interactions and other couplings. The renormalized single particle energies and 
the induced interactions are free of divergencies. The actual values of the renormalized quanti- 
ties depend on the special choice of the infinitesimal unitary transformation, i.e. on the cutoff 
A and on the special choice of the function u(x). On the other hand it is clear that measurable 
quantities like expectation values of observables should not depend on A or u(x). This is not a 
contradiction, since the unitary transformation applied to the Hamiltonian has to be applied to 
the observables as well. We have only calculated the effective Hamiltonian, thus we are not able 
to check this explicitely. But we have shown that the structure of the Hamiltonian does not de- 
pend on u(x). The renormalization of the single particle energies shows no essential dependence 
on u(x) if the cutoff A is sufficiently small. With some additional weak condition on u{x) (u(x) 
has to decay sufficiently fast), the induced interaction between two electrons forming a Cooper 
pair is attractive in the whole parameter space. 

Comparing our results for the induced electron-electron interaction with Frohlich's result 
( |1.3| ) shows that the effective interaction obtained by Frohlich has not only the problem that it 
contains a divergency, even the sign of the interaction is wrong in some part of the parameter 
space. In contrast the effective interaction obtained by Lenz and Wegner shows the correct be- 
haviour. Another difference is that in the perturbative expression of Frohlich the unrenormalized 
single particle energies enter. As a consequence, the renormalization of the electronic density 
of states, which has an important effect when one calculates T c , is not taken into account. As 
already mentioned in the introduction, a similar behaviour has been observed in the case of 
induced interactions for other problems as well [II]. The fact that the induced interaction in 
the general form ( |3.25 ) is always attractive can be traced back to the correct treatment of the 
different energy scales. Whenever large energy differences are treated first, the interaction (3.25) 
is attractive. If, as in the treatment by Frohlich || or by Bardeen and Pines Q all energy scales 
are treated at the same time, fk,-k,q,\ does not depend on k and q. In that case ( p. 25 ) yields a 
result similar to the one obtained by Frohlich. 

The discussion of the simple model in section 4 has shown, that the similarity renormalization 
scheme can be successfully applied to the electron-phonon problem. It yields the correct energy 
scale for single particle excitations and an effective interaction that yields reasonable values for 
the energy gap or for T c . But several problems have been left open. The reason was that the 
simple example was only presented to illustrate how the method works. Applying the similarity 
renormalization scheme to a more realistic model is possible and in that case a more detailed 
study should be done. For instance, the renormalization of the phonon frequencies has not 
been take into account. Furthermore, the renormalization procedure has been performed only 
to second order in the electron-phonon coupling. The similarity renormalization used here is a 
perturbative renormalization procedure. Higher corrections to the electron-electron interaction 
can be calculated systematically. The non-linear electron-phonon interaction in ( |3.I9| ) e.g. yields 
a contribution to the induced electron-electron interaction in fourth order. A detailed discussion 
of such higher terms may be important, since the attractive interaction is known to be a marginal 
relevant operator in a problem of interacting fermions. Furthermore it is possible to include 
the initial Coulomb repulsion in the Hamiltonian. It will be modified by the electron-phonon 
interaction as well, but it does not introduce a new energy scale to the problem. In this way 
it should be possible to obtain a more complete expression for the electron-electron interaction. 
We leave these problems for future work. 

Let us now compare the similarity renormalization scheme with Wegner's flow equations 
for Hamiltonians. A common feature of both approaches is that they work in a Hamiltonian 
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framework. As a consequence, we will never be able to obtain effective interactions containing 
retardation effects. The Hamiltonian is always a hermitian operator and induced interactions 
never depend on frequency. This is in contrast to methods using path integrals. In a path 
integral formulation of the present problem one is able to integrate out the phonons completely. 
The effective theory then contains a time-dependent or frequency-dependent electron-electron 
interaction. The interaction obtained in this way has nothing to do with the effective interaction 
in our approach. But this does not mean that one is not able to describe the decay of metastable 
states or the lifetime of some bound state in a Hamiltonian framework. In a recent paper on 
dissipative quantum systems this point has been discussed in great detail 13]. 

Both approaches, flow equations and the similarity renormalization scheme use a continuous 
unitary transformation to construct an effective Hamiltonian from a given initial Hamiltonian. 
This is typical for renormalization procedures. Already in the old formulation of renormalization 
Wilson [21] used a sequence of unitary transformations to calculate an effective Hamiltonian 
with a small cutoff. But the goal is different in both cases. In the flow equation approach, the 
Hamiltonian is transformed into a block-diagonal form. In the present case of the electron- 
phonon coupling the number of phonons is conserved in each block. But the Hamiltonian may 
still contain off-diagonal elements between two states with a large energy difference. In impurity 
models as in the Anderson impurity model [11|, or in dissipative quantum systems, where a 
small system is coupled to a bath [13], this does not happen, but there is no way to exclude 
such matrix elements in general. On the other hand, Wegner's flow equations can be used to 
diagonalize a given Hamiltonian approximately [10]. If in that case one performs the integration 
of the flow equations only to a finite value of the flow parameter, one obtains a Hamiltonian 
with a band-diagonal structure. The far off-diagonal matrix elements do not vanish, but they 
become exponentially small. 

When a Hamiltonian is transformed into a block-diagonal form, some interactions like the 
electron-phonon coupling in our case are eliminated completely. It is even possible to choose the 
continuous unitary transformation in such a way that some of the possible new induced inter- 
actions do not occur. This never happens in the similarity renormalization scheme. It has been 
designed to obtain an effective Hamiltonian with a band diagonal structure. As a consequence, 
the final Hamiltonian will contain any possible induced interaction. An example is the coupling 
c 9j a that occurs in our effective Hamiltonian. It has been avoided explicitely in the paper by 
Lenz and Wegner. On the other hand, all the induced couplings in the similarity renormalization 
scheme are free of divergencies. As a consequence, a systematic expansion of matrix elements 
in a coupling constant is always possible and well defined. The induced interactions one obtains 
in the flow equation approach are less singular than similar expression obtained perturbatively, 
but they may still contain some weak divergencies. An expansion in a small coupling constant 
is therefore not automatically well defined. But due to the simpler structure of the final Hamil- 
tonian, it is much easier to calculate expectation values of observables or dynamical properties 
in the flow equation approach. Therefore, depending on the problem one wants to treat, one or 
the other approach will have advantages. 
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Figure captions 



Fig. 1. The energy gap as a function of uj for 2#QiV(0)/u;o = 1/3 (a) and Ig^N (<S) / uq = 1.4 (b) 
and different values of T. In Fig. la the temperatures are T = (solid line), T = 0.0075u;o 
(dashed line), and T = O.Olwo (long dashed line). In Fig. 2a the temperatures are T = 
(solid line), T = 0.12cjo (dashed line), and T = 0.15u;o (long dashed line). The scaled curves 
(using Q4.8D) lie on top of the solid line. 



Fig. 2. The gap at to = as a function of temperature for 2gQN(0)/ujQ = 1/3 (a) and 2gQN(0)/ujQ 
1.4 (b). 



Fig. 3. T c as a function of the coupling 2gQN(0)/ujQ. The solid line shows the result from (4.6). The 



dashed line is the result obtained with the interaction of Lenz and Wegner, as explained in 



the text. The long dashed line shows the result of BCS-theory, obtained from (4.7), and 
the dashed-dotted line is a plot of the McMillan-Dynes formula. 



Fig. 1a 




Fig. 1b 



0.40 




0.020 



Fig. 2a 



0.015 



0.010 



0.005 



0.000 




0.000 



0.005 



0.010 



77(fl n 



